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ABSTRACT 



A computer program, with the code name TRUMP, was developed in 1965 
by the Lawrence Livermore Laboratory to solve a variety of transient 
and steady-state conduction heat flow problems in simple or complicated 
geometric configurations. The purpose of this investigation was to 
adapt "TRUMP" to the Naval Postgraduate School IBM/360 Model 67 com- 
puter system. Several heat conduction problems are solved with the 
adapted version of TRUMP and the results compare closely with the 
analytical solutions to these problems. The example problem inputs 
may be used as guides in preparing input data for future problems. 

Some suggestions are also given for further development of the program. 
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I. INTRODUCTION 



A computer program with the code name TRUMP, was developed in 1965 



solve a wide variety of non-linear, conduction heat transfer problems. 
TRUMP has been extensively used and continually improved since its 
initial development and has been successfully adapted to a number of 
computer systems including the CDC/6600 (LLL) , IBM/7094 (Bell Labs), 
UNIVAC/1108 (Battelle Northwest, Gulf General Atomics, NBS) , and the 
IBM/360 (B.F. Goodrich, NASA Lewis Research Center, Bendix Electro- 
dynamics, Cockerill-Ougee-Providence) , to name a few. The objective 
of this thesis project was to adapt TRUMP to the Naval Postgraduate 
School (NPS) IBM/360 Model 67 computer system. 

The B.F. Goodrich Corporation supplied a listing of a version of the 
TRUMP program that they had converted for use in their IBM/360 system. 

Cards were punched from this listing and a version of TRUMP was thus 

*■ 

secured . 

The B.F« Goodrich version was not complete. In addition to errors 
that resulted from punching 3,350 cards, the program also had an align- 
ment error in one subroutine of the program. Also, several options 
that are available in the complete LLL version of TRUMP were missing 
in the B.F 0 Goodrich version. 

With a great deal of assistance from the NPS Computer Facility staff, 
the majority of these problems have been corrected and a running version 
of TRUMP is now available at NPS. 

No attempt is made in this thesis to give a complete description of 
the TRUMP program. Rather, limitations of the NPS IBM/360 version as 



by Edwards 




of the Lawrence Livermore Laboratory (LLL) to 
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compared to the LLL CDC/6600 version are discussed. In addition, several 
example problems are solved and may be used as guides in preparing input 
data for future problems. For a complete description of the program, 
reference should be made to the TRUMP users manual 1 Reference 1 . 



t 
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II. DESCRIPTION OF THE PROGRAM 



A. GENERAL 

The program solves the general, non-linear, parabolic partial dif- 
ferential heat equation. Transient and steady-state heat conduction 
problems in multidimensions may be solved. Initial conditions, material 
properties, source and sink strength and boundary conditions may vary 
with time, spatial position and/or temperature. 

In order to solve problems with TRUMP, the user specifies input 
data in "BLOCK" form to 12 input data blocks. A description of each 
of these blocks and their use will be given in Section B 2. 

Following the definition of the problem a model must be constructed. 
Complex geometric shapes are divided into regions and if there is sym- 
metry, the symmetry planes are represented by insulated surfaces in 
order to simplify the model. Regions may have different materials with 
specified properties. 

The regions may be subdivided into volume elements or nodes having 
arbitrarily any size and shape. Nodes or volume elements must have a 
representative nodal point whose location may be anywhere in the node 
or on the surface of the node. However, in transient problems, for 
maximum accuracy, node shapes and nodal point locations are chosen so 
that the lines joining the nodal points of connected nodes are perpen- 
dicularly bisected by the connected area. In steady-state problems, 
due to the lack of heat capacity effects, the nodal points may be lo- 
cated anywhere in the node without loss of accuracy. A variety of 
boundary conditions may be given to the surface nodes. Constant temper- 
ature, variable temperature, constant heat flux, variable heat flux, 
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forced convection conditions, free convection conditions, radiation 
boundary conditions or combinations of the above may be treated. 

Instead of having very fine zoning of the regions, rather crude 
zoning will give the solution with less effort. Later, the solution 
may be obtained for a finer subdivision of the system. Also, with 
the crude zoning one can easily understand which regions need finer 
zoning. 

The accuracy of the given problem depends on how well the problem 
can be modeled, and the fineness of the spatial subdivisions and time 
steps used in the calculations. 

Allowable problem size limits for the NPS IBM/360 version of TRUMP 
are summarized in Table 1. 
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TABLE I 



Problem Size Limits 



Item 


Data Block 


Parameter 


Size 


Materials 


2 


• M 2 


15 


Reactants 


3 


m 3 


5 


Nodes 


4 


M 4 


355 


Internal thermal connections 


5 


m 5 


950 


External thermal connections 


6 


m 6 


60 


Boundary nodes 


7 


m 7 


20 


Heat generation tables 


8 


Mg 


5 


Initial conditions 


9 


m 4 


355 


Mass flow connections 


10 


M 10 


50 


Temp, versus time plot nodes 


11 


M 11 


10 


Temp, versus plot times 


- 


Mi 


1 


Remotely dependent properties 


12 


M 12 


75 


Table lengths 


2, 3, 6 , 7 , 8 , 10 


m q 
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In general, any consistent set of units may be used for the input 
data; output quantities are then given in the same unit system. 

Each material, system node, and boundary node can be identified 
by a non-zero integer number fitting within 5-column fields on the 
data cards. These identifications can be in any sequence. In order 
to reduce the number of input cards arithmetic sequencing is suggested 
so that the nodes having identical descriptions can be submitted as 
input on only one card. 

For non-symmetr ic , symmetric with respect to an axis, or to a 
center, the nodal volumes and the thermal connection areas can be cal- 
culated by a control value in the input data. 

Initial conditions may be specified in BLOCK 1 of the input data. 
These values will be assigned to all nodes described in BLOCK 4, unless 
other values are specified for individual nodes in BLOCK 9. 

Each internal thermal connection between nodes must be described 
in the input data by specifying the two node identification numbers, 
two connector lengths, and two interface dimensional factors. The 
selection of connector lengths depends on the location of the nodal 
point within each node. These locations may be arbitrary, but in 
transient problems should usually be at the geometric centers of the 
nodes. If the amount of the heat flow depends on surface temperatures, 
then they can be located on the surface. 

A system that exchanges heat with its surrounding must have external 
temperatures specified along its boundary. In the input data, boundary 
node identification numbers and their temperatures, which may be con- 
stant or a function of time, must be given. The term "boundary node" is 
identical in meaning to an external temperature. 
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Surface temperatures can be specified by using boundary nodes with 
specified temperatures and connecting the surface nodes to a boundary 
node with a very large convection coefficient. In order to reach a 
high degree of accuracy, surface nodes should be selected as either 
zero-volume nodes, or very thin nodes. When a large number of surface 
nodes are connected to the same boundary node and the heat transfer 
coefficients are not tabulated, the number of external connections can 
be reduced by replacing them with connections to a single zero-volume 
node that is connected to the boundary node with a very large convection 
coefficient . 

In order to generate steady-state solutions, TRUMP must first go 
through a set of transient calculations. The steady-state solution is 
obtained by allowing enough problem time to take the transient to its 
conclusion. However, if only steady-state results are desired, the 
control on the average temperature change during a time step (TVARY) 
can be made very large. This allows the program to neglect transient 
stability limitations on its choice of time steps resulting in a mini- 
mum expenditure of computer time on transient calculations. 

B. USE OF COMPUTER 

1 . General 

The data deck must include a problem name card, any number of 
block number cards X'/ith their input cards, and a data end card. The 
first card of each data deck must be a problem name card. The symbol 
must be in column 1, any identification of the problem in columns 2 
through 71 and 73 through 80. 
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In the program there are 12 allowed input data blocks, with 
block numbers from 1 through 12 „ Each block must begin with a block 
number card, and except for BLOCKS 1 and 11, must end with a blank 

card. The block number card must begin with the word "BLOCK" in 

columns 1 through 5, and the block number, from 01 through 12, in 
columns 6 and 7. 

Data BLOCK 1 must follow the problem name card; all other data 
blocks may be placed in any order. 

The last card of each data deck must be a data end card, with 

the word "ENDED" in columns 1 through 5, and "-1" in columns 6 and 7. 

Additional cards containing comments without "*" in column 1 
may precede the name card and will appear on the printout. Any desired 
block description in columns 9 through 80 may be made on the block 
number card. This will also appear on the printout. 

2 . Block Item Descriptions 

Basically, one should refer to the TRUMP users manual £ Reference 1 
for format details and the original descriptions of the input values. 

In this section the main difference between the IBM/360, Model 67 version 
of TRUMP and the LLL, CDC/6600 version are discussed. 

In BLOCK 1; 

MCYC Maximum allowed number of time steps. If zero, neg- 
ative, or unspecified will be 30,000. 

MSEC Maximum allowed machine time, in seconds. If zero, 
negative, or unspecified will be 30,000. 

NPUNCH Indicates that the final values of temperatures of nodes, 
weight fractions of reactants, and constant heat gener- 
ation rates of nodes will be punched. Since punched 
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output is not obtained in this version of TRUMP, the 
value of NPUNCH should not be used. 

IRITE Number of time steps between data output on unit ITAPE 
of problem time and table of node numbers and node 
temperatures, in addition to output of first, second, 
and last time steps. Only used in TRUMP/6600/DS version. 
This control is not in use in this version. 

ITAPE Output unit on which output data will be written as 
described in IRITE. This input value is not in use 
for IBM/360 Model 67 version of TRUMP. 

TIMEP Problem time interval between data output, in addition 
to output on first, second, and last time steps. This 
input control value is also not in use for this version 
of TRUMP. 

In addition to the above differences, the data transfer option 
is not in use in this version of TRUMP. 

3 . Numerical Output. 

The amount of output data is controlled by the value of KDATA. 
However, for every problem, all input data and several other quantities 
such as node volumes and connection areas are written out. The results 
of the first, second, and last time steps are always written out. In 
addition, output data are written out at problem time step intervals 
of IPRINT. For the notations used on the printouts one should refer 

-- -j 

to TRUMP users manual [Reference 1 . 

4 . Evaluation of Results 

Any desired end to the problem may be made with the specified 
input values in BLOCK 1. These may be maximum problem time,TIMAX, 
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maximum temperature, TMAX, minimum temperature, TMIN, maximum number 
of time steps, MCYC, or maximum number of seconds of machine time, 
MSEC. On the printouts, the value of the problem-end sentinel, KWIT, 
will appear. Table II shows the values of KWIT, and the descriptions 
of these values . 
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KWIT 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 



TABLE II 



Values of KWIT 



Cause of Problem Ending 

The problem time, SUMTIM, reached TIMAX 

A temperature exceeding TMAX 

A temperature less than TMIN 

Steady-state 

Data 

Relates to the tape usage 

The number of time steps reached MCYC 

Machine time reached MSEC 

BLOCK 2 or BLOCK 4 missing 

Convergence failure 

Items exceeded the size limits 

Table lengths exceeded the size 
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C. DESCRIPTIONS OF SUBROUTINES 



THERM 

TALLY 

CHEM 

SPECK 



FINK 

GEN 

SURE 



FLOW 

SEEK1 

SEEK2 



Subroutine for material properties and node descrip- 
tions. BLOCKS 2, 4, and 12 are included. 

Subroutine for initialization, totaling, and checking 
results, converting regular to special nodes, testing 
for ending the problem. BLOCKS 1 and 9 are included. 
Subroutine for chemical reaction. BLOCK 3 is included. 
Subroutine for heat flow between special nodes and 
other nodes or boundary nodes . Iterates to solve 
set of implicit difference equations when special 
nodes connected to each other. 

Subroutine for heat flow between nodes by conduction, 
convection, and radiation. Treats all nodes as regu- 
lar nodes. BLOCK 5 is included. 

Subroutine for internal heat generation. BLOCK 8 
is included. 

Subroutine for heat flow between surface nodes, bound- 
ary nodes by radiation, free and forced convection 
treats all surface nodes as regular nodes. BLOCKS 6 
and 7 are included. 

Subroutine for mass flow between nodes. Treats all 
nodes as regular nodes. BLOCK 10 is included. 
Subroutine for finding an identification number of 
a material, reactant, node or boundary node in a 
block item list. 

Subroutine for finding the identification numbers 
of a pair of nodes in a block item list. 
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CLOCK, CLOCKl Subroutines for initialization before first time 



PLOT 



PATCH 



REFER 



SRCON 



step . 

Subroutine for making CRT plots and tables of temp- 
erature versus node location, temperature versus 
time, and node location versus time. BLOCK 11 is 
included . 

Subroutine for converting a number with a decimal 
point to floating number, otherwise substituting 
a specified floating point number. 

Subroutine for finding array subscripts of materials, 
reactants, nodes or boundary nodes referred to by 
identification number in another data block, and 
writing out diagnostic statements whenever one can 
not be found. 

Subroutine which is related to PATCH. Converts a 
number having decimal point to a floating point 
number . 
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III. EXAMPLE PROBLEMS 



A. SLAB -MELT PROBLEM 

A semi-infinite slab (-©o<^x <(oo , - oo <y < CxO , 0 C z < C o ° ) > 
has thermal conductivity 1.0 cal/sec-°C, density 1.0 g/cm heat capacity 
1.0 cal/sec-°C, a melting point of 50.0°C, and a latent heat of fusion 
of 25.0 cal/g. The initial temperature is 100. 0°C, and for times t^>0, 
the surface, z = 0, is maintained at 0.0°C. TRUMP is used to calculate 
the temperature distribution in the slab, and the motion of the solid- 
liquid interface, Zint> for times up to 200 sec. The results are com- 
pared with the exact solution. 

1 . Analytic Solution 

From Reference 3 analytic results are: 

Z 

= 112.4 erf ( ) for 0 <^T < 50°C 

2 ft 

Z 

T 2 = 9.9 + 90.1 erf (• — ) for 50 <^T ^ 100°C 

2 ft 

Zfnt = 0.835* IT for T = 50°C 

2 . Trump Solution 

The input data for this problem are shown in Figure 1. In order 
to have output at every 10 cycles, in addition to output on the first, 
second, and last time steps, IPR.INT is selected as 10 in BLOCK 1. For 
temperature accuracy, TVARY is 1.0°C and the problem is ended at 200 sec 

by setting TIMAX equal to 200.0. TONE is 100.0 since the initial tem- 
perature of the slab is 100. 0°C. However, since the zero-volume sur- 
face node, node -1, is at 0.0°C for t 0 its temperature is specified 
in BLOCK 9 by setting the value of TT equal to 0.0. In BLOCK 2 y the 
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material name, AMAT, and material identification number are arbitrarily 
selected at SMELT and 50, respectively. Other slab thermal properties 
such as density, heat capacity, thermal conductivity, temperature at 
which the latent heat is released and the latent heat of fusion are also 
specified in BLOCK 2. In order to give node descriptions in BLOCK 4 
the slab is divided into nodes as shown below. 



C\j 



1cm 



2cm 



12 



v- 

-+Z 
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--N-- 



-V-- 



4cm 



k- 



8cm 



31 



32 



-v 



V 



— *> 




• 


! 

• 


35 


36 



A total field depth of 100.0 cm is chosen so that, at 200 sec no 
appreciable temperature change will take place at the maximum value 
of Z. In this figure, the node identification numbers, and the node 
thicknesses are shown. The nodal points of all nodes are at the node 
centers. In BLOCK 5, internal connection node numbers, lengths of 
heat conduction paths from the nodal points to the connected interfaces, 
and the area of the connected interfaces are defined. In BLOCK 6, 
a zero-volume surface node, node -1, is connected to a boundary node, 
node 2001, with a heat transfer coefficient of 1.0 x 10^cal/sec-cnk-°C , 
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and the external surface area of the surface node is selected as unity. 



In BLOCK 
Computer 



7, a constant boundary node temperature is defined as 0.0°C. 
running time for this problem was 65.52 sec. 
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FIG. 1 INPUT DATA FOR SLAB MELT 



100 . 
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FIG. 2 TEMPERATURE DISTRIBUTION IN THE SLAB-MELT PROBLEM 
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DISTANCE (cm) 



TABLE III 



SOLID-LIQUID INTERFACE IN THE SLAB MELT 
Time (sec) z int (Analytic) z int (Trump) 



3.56 


1.57 


1.37 


4.19 


1.71 


1.77 


4.77 


1.82 


2.00 


5.92 


2.03 


2.00 


8.13 


2.38 


2.04 


9.67 


2.60 


2.46 


10.81 


2.74 


2.86 


12.08 


2.90 


3.00 


15.90 


3.33 


3.02 


18.43 


3.58 


3.44 


20.13 


3.75 


3.84 


22.21 


3.94 


4.00 


28 . 24 


4.44 


4 . 16 


30.91 


4.64 


4.56 


33.10 


4.80 


4.97 


38.20 


5.16 


5.00 


44.21 


5.55 


5.37 


47.20 


5.74 


5.77 


50.42 


5.93 


6.00 


60.16 


6.48 


6.23 


64.03 


6.68 


6 . 63 


67.34 


6.85 


7.00 


78.78 


7 .41 


7.13 


83.62 


7.64 


7.54 


87.67 


7.82 


7.94 


106.12 


8.60 


8.46 


110.92 


8.79 


8.87 


120.92 


9.18 


9.00 


131.54 


9.58 


9.42 


137.10 


9.78 


9.82 


146.27 


10.10 


10.00 


160.05 


10.56 


10.39 


166.34 


10.77 


10.79 


175.46 


11.06 


11.00 


191.54 
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B. FINITE-SLAB PROBLEM (suggested by Professor Paul Pucci) 

A large, 1.0 inch thick steel plate, initially at 80.0°F, is 
suddenly submerged in a fluid at temperature 680. 0°F. TRUMP is used 

to calculate the transient temperature distribution in the plate assuming 
a constant surface heat transfer coefficient of 1248.0 BTU/hr-f t 2 -°F. 

The thermophysical properties of the steel were taken as; thermal con- 
ductivity 26.0 BTU/hr-f t-°F, density 487.0 lbm/ft'f and specific heat 

0.133 BTU/lbm-°F. 

1 . Analytic Solution 



(x,t) 



T - T 



oO 



= 2 



Ti - T_ 



SinAj^L 



-(A n L) 2 Fo 



, A R L + SinA n LCos A n L 



Cos(A n L*) 

Ll 



hL 

Applying the Biot number, Bi = — ~ 2.0 to the eigen value 

A n L 

problem of Cot(A n L) = j the two eigen values AiL = 1.0769 and 

Bi 

A^L = 3.6436 can be found. After substituting the appropriate value of 

m oC t 

the Fourier number, Fo = 2 = 230. 4t in the above series solution and 

L 



considering only the first two terms, the result becomes 
= 680 - 1200 0 .58923153 e 267 ' 2t Cos (1.0769*) - 0.11836085e 



-3058. 75t „ 

Cos (3 . 6436jp) 

w 



2 . Trump Solution 

The input data for this problem are shown in Figure 3. For unit 
consistency, lengths in inches are converted to feet and applied to the 
problem. Since the plate has symmetry with respect to the mid-plane, 
this plane is taken as an insulated plane in the TRUMP solution. So, 
in BLOCK 1 symmetry type indicator, KD, is applied as 1 for the non- 
symraetric case, and the total thickness of the steel plate is considered 
as 0.5 inches. Temperature scale indicator, KT, is selected as 3 
because of the Fahrenheit scale. For temperature accuracy, the desired 
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maximum temperature change in each time step, TV ARY, is selected as 
1.0°F. Maximum problem time, TIMAX, is given as 1.0 hr. TONE is set 
equal to 80.0 making the initial temperature of all nodes equal to 80.0°F. 

In BLOCK 2, material name and material number are given as, FSLAB and 
50, respectively. Also, the steel properties such as density, specific 
heat, thermal conductivity are defined. In BLOCK 4 node numbers, mate- 
rials of nodes and their dimensions are given with respect to the equally 
spaced 10 nodes for the 0.5 inches of thickness of the slab. In BLOCK 6, 
the zero-volume surface node, node -1, is connected to a boundary node, 
node 2001, with a surface heat transfer coefficient of 1248.0 BTU/hr-ft^-°F. 
In BLOCK 7, a boundary node temperature of 680. 0°F is given. The com- 
puter running time for this problem was 24.46 sec. 

In order to see the effect of finer zoning on the results, the 
same problem was also solved using 20 and 40 nodes. However, the results 
were all within 1°F accuracy and not much improvement was obtained with 
the finer grid sizes. Computer running time, on the other hand, was 
increased from 24.46 sec. to 53.76 sec. and to 172.64 sec. for 10, 20, 
and 40 nodes, respectively. 
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FIG. 3 INPUT DATA FOR FINITE SLAB 



TEMPERATURE (°F) 




NODE NUMBERS 

FIG. 4 TEMPERATURE DISTRIBUTION IN THE FINITE SLAB PROBLEM 
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TEMPERATURE ( °F) 




100 



_l_ TRUMP Solution. 



0 . 1 . 2 , 3 . 4 . 

T IME (sec) 



FIG. 5 SURFACE TEMPERATURE HISTORY OF FINITE SLAB PROBLEM 
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C. SEMI- INFINITE SLAB PROBLEM 



A semi-infinite slab is initially at a uniform temperature of 100. 0°C. 
At time t = 0, the surface temperature is suddenly lowered and main- 
tained at a temperature of 0.0°C. The slab has a thermal conductivity 



heat flow from the slab are calculated as a function of time for times 
up to 200 sec using TRUMP, and the results are compared with the ana- 
lytical solutions. 

1 . Analytical Solution 



2 . Trump Solution 

The input data form is shown in Figure 6. In BLOCK 1, problem 
controls such as problem time, TIMAX = 200.00 sec, maximum temperature 



all nodes, TONE = 100. 0°C are given. In BLOCKS 2, 4, 5, 6, 7 and 9, the 
thermal properties, zoning procedure, identification numbers, internal 
and external connection descriptions, boundary node temperatures, and 
initial conditions are defined identically as in the slab-melt problem. 
Computer running time for this problem was 27.05 sec. 



of 1.0 cal/sec cm-°C, density of 1.0 g/cm^, and a heat capacity of 
1.0 cal/sec-°C. The temperature distibution in the slab, and the total 



T (x , t ) = 100 erf 




4t 




change in each time step, TVARY = 1.0°C, and the initial temperature of 
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FIG. 6 INPUT DATA FOR SEMI-INFINITE SLAB 
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FIG. 7 TEMPERATURE DISTRIBUTION IN THE SEMI- INFINITE SLAB PROBLEM 
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TABLE IV 



HEAT FLOW FROM SEMI- INFINITE SLAB 



Time (sec) 


Qtot (Trump) 


Qtot (Analytic) 


1.26 


-119.42 


-126.56 


1.81 


-146.12 


-151.86 


2.65 


-178.93 


-183.56 


3.85 


-217.69 


-221.43 


5.63 


-264.69 


- 267.69 


8.21 


-320.90 


-323.39 


12.06 


-389.63 


-391.79 


17.70 


-472.82 


- 474.75 


25.97 


- 573.36 


-575.00 


38.02 


-694.20 


- 695.73 


55.75 


-841.12 


-842.50 


81.76 


-1019.00 


-1020.29 


119.93 


-1234.40 


-1235.71 


175.44 


-1493.20 


-1494.59 


200.00 


-1594.30 


-1595.77 
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TRUMP OUTPUT DATA DATA DECK 
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D . CYLINDER PROBLEM 



An infinitely long cylinder of unit radius, with a unit surface 
temperature, is embedded in a medium with a zero initial temperature. 
TRUMP is used to calculate the temperature distribution in the medium. 
All thermal properties for the medium have unit values. Results are 
compared with an analytical solution. 

1 . Trump Solution 

The input data are shown in Figure 8. Because of the geomet- 
rical shape of cylinder, the symmetry type indicator, KD, in BLOCK 1 
is selected as 2. Maximum problem time, TIMAX, is made 50, and the 
medium initial temperature, TONE, is set to 0.0. Also, for temperature 
accuracy, TVARY, is selected as 0.05. In BLOCK 2, unit properties of 
the medium are defined. Nodal point locations are chosen so as to be 
able to compare the TRUMP solution with the given table in Reference 5 
rather easily. The zoning of the medium is shown in the figure below. 
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In BLOCK 4, the arithmetic mean average radii, DRAD, are used to obtain 
correct nodal volumes for cylindrical shaped nodes. Actual node inter- 
faces are used as interface radii, DRAD, in BLOCK 5, to obtain the areas 
of the connected interfaces. In BLOCK 6, zero-volume surface node, 
node -1, is connected to the boundary node with a heat transfer coeffi- 

O 

cient of 1.0 x 10°, and geometric factors for external surface node 
area are defined. Computer running time for this problem was 15.46 sec. 
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INPUT DATA FOR CYLINDER 
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ENDED-1 LAST CARD OF DATA DECK 



TABLE V 



TEMPERATURE DISTRIBUTION FOR CYLINDER (TIME = 50) 



R T (Trump) T (Analytic) 



1.0 


1.0 


1.0 


1.1 


0.961 


0.963 


1.3 


0.897 


0.898 


1.5 


0.842 


0.843 


1.7 


0.794 


0.794 


1.9 


0.751 


0.751 


3.0 


0.574 


0.574 


5.0 


0.382 


0.381 


7.0 


0.261 


0.260 


9.0 


0.179 


0.177 


10.0 


0.146 


0.146 


20.0 


0.016 


0.016 


30.0 


0.001 


0.001 
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TRUMP OUTPUT DATA DATA PECK 
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E . SPHERE PROBLEM 



A sphere with a unit radius and a unit surface temperature is 
embedded in a medium with a zero initial temperature. The temperature 
distribution in the medium is calculated using TRUMP assuming all the 
thermal properties have unit values. Results are compared with the ana- 
lytical solution. 

1 . Analytic Solution 

r-1 

T(r,t) = erfc(-— ) 

2 . Trump Solution 

The input data from the cylinder problem can be used to solve 
this problem if the symmetry type indicator, KD, in BLOCK 1 is changed 
from 2 to 3 to correspond to the centrisymmetric geometry. All other 
input conditions are identical with the cylinder problem. The computer 
time used for this problem was 11.61 sec. 
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TABLE VI 



TEMPERATURE DISTRIBUTION FOR SPHERE (TIME = 50) 



R T (Trump) T (Analytic) 



1.0 


1.0 


1.00 


1.1 


0,89 


0.90 


1.3 


0.75 


0.75 


1.5 


0.64 


0.64 


1.7 


0.56 


0.56 


1.9 


0.49 


0.49 


3.0 


0.28 


0.28 


5.0 


0.14 


0.14 


7.0 


0.08 


0.08 


9.0 


0.05 


0.05 


10.0 


0.04 


0.04 


20.0 


0.00 


0.00 


30.0 


0.00 


0.00 
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F. TOO DIMENSIONAL PLATE PROBLEM 



A rectangular plate, 40 x 20 cm, has a thermal conductivity of 
400 cal/cm-sec-°C, a density of 1.0 g/cm^, and a heat capacity of 
1.0 cal/g-°C and is initially at 0.0°C. Three sides of the plate are 
maintained at 0.0°C. The upper side of the plate has a temperature 
distribution that varies sinusoidally with x, starting at 0°C, reaching 
a maximum of 100°C at the midpoint and decreasing to 0°C at the outer 
edge (see figure below). The steady-state temperature distribution is 
calculated in the plate using TRUMP. Results are compared with the 
analytical solution. 




1 . Analytic Solution 



T(x,y) = 100 



-Ry 

Si nh ( ~~) 
S inh ( ) 



Sin(-'^- ) 
40 



where 0 ^x ^ 40 and 0 ^y ^ 20 
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2 . Trump Solution 



The input data are shown in Figure 9 . Because the boundary 
conditions of the plate and the plate itself are symmetric with respect 
to the plane at x = 20 cm, this symmetry plane is assumed to be an in- 
sulated surface. The zoning of the plate is shown in Figure 10. In 
BLOCK 1 of the input data, IPRINT, is selected as 9999, so that print- 
outs are obtained only for the first, second, and last time steps. To 
obtain a rapid conclusion to the transient problem, TVARY, is selected 
as 100. 0°C. In BLOCK 4, instead of using a very large number of surface 
nodes for sides x = 40 cm and y = 0, surface node -1 and surface node -4 
are used, respectively, in order to reduce the number of the input cards, 
and in BLOCK 6 they are connected to the boundary nodes 2001 and 2004. 

In BLOCK 7, the boundary node temperatures, and in BLOCK 9, the initial 
temperatures of the surface xxodes are giveix. The computer runnixxg time 
for this problem was 28.71 sec. 
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TWO-DIMENSIONAL SLAB PROBLEM. NOVEMBER 1971 
CK 1 PROBLEM LIMITS AND CONTROLS. 
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FIG. 9 INPUT DATA FOR TWO- DIMENSIONAL -PLATE 
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FIG. 10 ZONING OF THE TWO DIMENSIONAL PLATE PROBLEM 
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TABLE VII 



TEMPERATURE DISTRIBUTION FOR TWO-DIMENSIONAL PLATE IN STEADY-STATE 



Node Number 


T (Trump) 


T (Analytic) 


1 


0.26 


0.27 


5 


2.17 


2.23 


10 


3.34 


3.41 


11 


0.79 


0.81 


15 


6.57 


6.71 


20 


10.10 


10.30 


21 


1.35 


1.37 


25 


11.14 


11.37 


30 


17.13 


17.45 


31 


1.93 


1.97 


35 


16.02 


16.31 


40 


24.63 


25.03 


41 


2.58 


2.62 


45 


21.32 


21.65 


50 


32.77 


33.24 


51 


3.28 


3.33 


55 


27.19 


27.53 


60 


41.77 


42.26 


61 


4.08 


4.12 


65 


33.75 


34.09 


70 


51.84 


52.33 


71 


4.97 


5.01 


75 


41.18 


41.49 


80 


63.24 


63.69 


81 


6.00 


6.03 


85 


49.65 


49.92 


90 


76.23 


76.62 


91 


7.17 


7.20 


95 


59.36 


59.58 


100 


91.12 


91.45 
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IV. CONCLUSIONS 



The adaptation of the computer program "TRUMP" to the Naval Post- 
graduate School IBM/360 Model 67 has been achieved, with success. 
Analytical temperature distributions of various transient and steady- 
state heat conduction problems compared very well with the results 
obtained from TRUMP. The complete adaptation has not yet been estab- 
lished, however. At present, punched output and plotted graphs have 
not been obtained. The computer solutions of the example problems did 
not require the use of subroutines CHEM, GEN, and FLOW. Thus, the 
adaptation of these subroutines is yet to be evaluated. Since, the 
accuracy of the solutions are dependent on the zoning, this effort can 
be reduced if an existing program "FED: A computer program to generate 
geometric input for the heat - trans f er code TRUMP" by Dale A. Schauer 
Reference 4 is adapted to the IBM/360 Model 67 version of TRUMP. 

The example problems presented herein will provide a useful reference 
for preparation of input data for future problems. 
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